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The ratio of the Zeeman splitting to the cyclotron energy (M = AEz/hiVc) for hole-like carriers 
in bismuth has been quantified with a great precision by many experiments performed during the 
past five decades. It exceeds 2 when the magnetic field is along the trigonal axis and vanishes in 
the perpendicular configuration. Theoretically, however, M is expected to be isotropic and equal 
to unity in a two-band Dirac model. We argue that a solution to this half-a-century-old puzzle can 
be found by extending the k ■ p theory to multiple bands. Our model not only gives a quantitative 
account of magnitude and anisotropy of M for hole-like carriers in bismuth, but also explains its 
contrasting evolution with antimony doping pressure, both probed by new experiments reported 
here. The present results have important implications for the magnitude and anisotropy of M in 
other systems with strong spin-orbit coupling. 


Spin-orbit interaction (SOI) is a central issue in con¬ 
temporary solid state physics. It is an automatic con¬ 
sequence of the Dirac theory, and well understood for a 
lonely electron in presence of a single atomic potential. 
In crystals, however, a diversity of SOI effects arise due 
to variety of the crystalline potentials. The effect also 
strongly depends on the momentum since the magnitude 
of SOI is set by (7i/4m^c^)cr • VY(r) x p. (Here V{r) 
is the crystalline potential, p is a momentum, and cr is 
the Pauli spin matrix.) It is thus basically difficult to 
study the effect of SOI for various materials and various 
/c-points of the Brillouin zone by a universal approach. 

One signature of crystalline SOI is its impact on the 
ratio of the Zeeman splitting AEz to the cyclotron en¬ 
ergy fiuic, dubbed M = AEz/huic- The crystalline SOI 
appears as an antisymmetric part with respect to commu¬ 
tation of momentum operators in the Hamiltonian under 
a magnetic field [ll-@|. Its eigenvalue corresponds to the 
effective Zeeman energy with an anisotropic effective g- 
factor g, while the eigenvalue of the symmetric part cor¬ 
responds to the cyclotron energy with an anisotropic ef¬ 
fective mass. Therefore, the impact of the crystalline SOI 
can be characterized by the relative energy scale of the 
crystalline SOI to the kinetic energy, i.e., M = AEz/Tt-oJc- 
This ratio can be accurately determined by experiment in 
those cases in which quantum oscillations simultaneously 
detect successive Landau levels as well as same-index sub- 
levels with opposite spins [Ml- The energy levels for 
different cases and their corresponding M are illustrated 
in Fig. [1] and the experimental values of M in three 
different systems are listed in Table [H As seen in this 
Table, when SOI is weak (i.e., in the case of graphite), 
M is much smaller than unity. One can show that M is 
exactly equal to unity (and so never exceeds unity) for 
any direction of magnetic field when a large SOI strongly 
couples two bands based on the two-band model, which 
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FIG. 1. Energy levels under a magnetic field for different 
ratio of the Zeeman splitting to the cyclotron energy, M = 
AEzIhujc- 

is equivalent to the Dirac Hamiltonian S Si- This 
agrees quite well with the experiments on the L-point of 
bismuth. 

However, in the case of holes at T-point, there are two 
puzzling features: (i) M is extremely anisotropic (~ 0 for 
B T trigonal axis); (ii) It largely exceeds unity in one con¬ 
figuration (= 2.12 for B II trigonal axis) as shown later. 
These two puzzles emerged from numerous experiments 
starting half a century ago [see supplementary materials 
(SM)] in absence of any satisfactory explanations [ffl. 
In other words, the two-band approach completely fails 
to give its adequate value. 

In this Letter, we show that a satisfactory solution to 
these longstanding puzzles can be found by going beyond 
the two-band Dirac model. Furthermore, we present new 
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TABLE I. Ratio of the Zeeman splitting to the cyclotron energy for some compounds. The cyclotron mass UM, [il,[a, 
effective g-factor , and atomic SOI [J are also listed. For graphite, only values with the magnetic field parallel to 

the c-axis are listed. The range of values for Bi expresses their anisotropy. 



Graphite 

InSb 

Bi 


position 

K (ele.) 

H (hole) 

r 

L (ele.) 

T (hole) 

mc/m 

0.038 

0.057 

0.014 

0.0019-0.027 

0.068-0.22 

9 

2.5 

2.5 

52 

74-1060 

0.79-63 

M = mcg/2m 

0.048 

0.073 

0.36 

0.9-1.0 

0.0-2.12 

atomic SO (eV) 

0.005 

0.005 

0.27, 0.68 

1.8 

1.8 


experimental results on the evolution of M for holes in 
bismuth with pressure and antimony doping. Rather 
counterintuitively, these two alternative ways to reduce 
carrier concentration shift the magnitude of M in op¬ 
posite directions. We show this can be quantitatively 
explained by the present approach. 


Numerous experiments have quantified the magnitude 
of the g-factor (and therefore M) for both electrons and 
holes in bismuth [Ml. The large atomic SOI of bis¬ 
muth (~ 1.8 eV) overwhelmingly dominates all other 
energy scales including the Fermi energy (~ 28 meV) 
and the band gap (~ 15 meV). Therefore, the crystalline 
SOI dramatically affects the electronic structure and this 
leads to a very complex hierarchy between Tuoc and Ai5z 
for different carriers and different orientations of mag¬ 
netic field. In the case of electrons at the L-point of 
the Brillouin zone, M ~ 1 with little dependence on the 
field orientation. In other words, in spite of the extreme 
anisotropy of both the cyclotron mass rric and g ll| , huJc 
and Aiilz remain almost equal to each other. This prop¬ 
erty of the L electrons is quite well-understood based on 
the two-band model and were employed to give 

a quantitative account of the complex Landau spectrum 
of L electrons 11, 1^ . The puzzling features experimen¬ 
tally observed on holes at the T-point are shown in Fig. [2] 
(a), which shows the angular dependence of two Nernst 
peaks corresponding to energy levels n = 2^, where ± 
indicates the degree of freedom of the Kramers doublet 
The angular dependence of M can be deduced by 
plotting — B^~), which is shown in Fig. [2](b), 

where is the magnetic field of the Nernst peak for 
n = 2^ and Ftri is the oscillation frequency for i? || trigo¬ 
nal axis. For B || trigonal, this value is exactly the same 
as M, and so we obtain M = 2.12 (cf. Fig. [T|). For B T 
trigonal, on the other hand, F~^{B‘^^ — B‘^~) becomes 
almost zero indicating M ~ 0. 


Next we give a general theory on the ratio M. We start 
from one-electron Hamiltonian in presence of strong SOI 
and apply k ■ p theory to a multiband (n-band) system 
under a magnetic field taking into account the SOI in a 
fully relativistic (non-perturbative) way. The obtained 
Hamiltonian is written in terms of 2n x 2n matrix for 
an n-band system. (See SM for details.) In order to 
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FIG. 2. Angular dependence of a Kramers doublet in bis¬ 
muth. (a) Magnetic field at which n = 2'^,2~ Landau levels 
are evacuated as a function of the orientation of the mag¬ 
netic field. 9 is the angle between the magnetic field and the 
trigonal axis. Symbols represent the field at which peaks are 
resolved in the Nernst quantum oscillation data (Ref. ill), 
(b) Angular dependence of (B^+ — ) from panel (a). 


study the electromagnetic properties of electrons belong¬ 
ing to a particular band, we decouple 2x2 Hamiltonian, 
J^=o, from the other bands by using the Lowdin’s uni¬ 
tary transformation [l8l |. The cyclotron energy is the 
eigenvalue of the symmetric part of J^=q. A straight¬ 
forward calculation yields hiOc = {eTiB /c)y^det 
where a is the inverse mass tensor given by 
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We have extended the notations of Cohen-Blount as = 
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and Un = Vq^^, where v^j is the interband matrix 
element of the velocity operator between i-th band with 
spin ^ and j-th band with spin rj. Ei is given as the 
eigenvalue of the antisymmetric part of For B \\ i, 

it is obtained as Ej. = where g = 2m^/Gii, 


Gii — 4 


X ' tyj X Un 

^ Eq - En 




I ^ ^ •'n ~E ^ 


Eq — En 


and /tb = eTi/2mc. Then, M for i? || z is 


( 2 ) 


/ G,, 

TiUJc y Otxx^yy ^“x 


( 3 ) 


The obtained results 0 -® are all gauge independent. 
It should be stressed that the denominator Eq — En of 
a and G are approximately canceled each other in M. 
Therefore, M is very sensitive to the symmetric proper¬ 
ties of the interband matrix elements and insensitive to 
the energy differences. If only two bands (n = 0,1) are 
taken into account, Gzz = ctxxCtyy — c?xyi then M 
would be exactly unity, consistent with previous results 
[3, B 0|- If one take into account more than two bands, 
one can see that the form of Eqs. 0-0 already tells us 
that M can be larger than unity. 

Various theories of g for multiband systems have been 
studied so far. The simple formula available at the 
present moment is valid only in the seiniclassical limit 
[l9j, whereas the formulae for Bloch bands based on the 
quantum treatment are too complex to compute g for 
various systems BUS [llj. The present quantum for¬ 
mulae for hu)c^ g, and M are general, rigorous within 
k ■ p theory, and yet easy to handle. It is the advantage 
of these formulae that these values can be automatically 
obtained from the interband matrix elements and the en¬ 
ergy differences, which can be directly computed by the 
band calculations as shown later. 

Here we give a concrete example of the above general 
theory. We adopt obtained formulae to the T-point holes 
in bismuth by taking into account the symmetry at the 
T-point. But it must be noted that the following argu¬ 
ments are also valid for group V semimetals (Sb, As), 
IV-VI narrow gap semiconductors (PbTe, PbSe, SnTe, 
etc.), and the topological insulator Bi 2 Se 3 , since the k- 
points where their carrier locate have the same symmetry 
as the T-point of bismuth 22, 2^. The hole band at the 
T-point has the symmetry of 2J, [2^. The symme¬ 
tries of the other bands at the T points are shown in Fig. 
m (e) and Fig. 1 in SM. (The group theoretical notation 
is that of Ref. [1^.) From the selection rules, the ma¬ 
trix elements between the T^ band and the others are 
finite only for t^n^ = (T 4 ^(l)|t;|Tg+(n)) = (-a„,ia„, 0 ), 
Un'’ = {T^^{1 )\v\GTq (n)) = (-a„,-ia„, 0 ), and Un^'’ = 
(T 4 ^(l)|t;|C'T^(n)) = (0,0, 6 „), where C is the prod¬ 
uct of space inversion and time reversal operators and 


a„, bn are complex numbers [1^. The x, y, and z- 
directions are taken along the binary, bisectrix, and trig¬ 
onal axes, respectively. Since g is given by the outer 
products of tn and [Fq. ([2])], it is clear from tn'’ 
and that Gzz is the only non-vanishing term. 

Therefore, the angular dependence of ATz is simply 
AEz{0) = 2m/iBT| cos0|i/G^, where 9 quantifies the 
tilt angle of the magnetic field off the trigonal axis. This 
leads to 5 = 0 for B T trigonal axis, providing a solution 
to the first experimental puzzle. 

The ratio M for the T-point holes is obtained as 





Eq — En 
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( 4 ) 


This form makes it easy to understand how M can ex¬ 
ceed unity. For example, in the case with three bands 
n = 0,1,2, whose energies are E 2 < Eq < Ei, it is 
easy to show that the ratio is always greater than one 
no matter how large the energy differences are [Fq. (32) 
in SM]. Accordingly, one expects that the interband con¬ 
tributions from the lower (higher) energy bands increase 
(decrease) the ratio, providing a possible solution to the 
second experimental puzzle. 

In order to know the definite value of M, we need to 
evaluate the interband matrix elements tm Un and the 
energy differences Eq — En based on the band calcula¬ 
tions. Here we evaluate them based on the multiband 
k ■ p Hamiltonian (eight-band model) derived from the 
well-known tight-binding band calculation of Bi by Liu 
and Allen, which is in quantitative agreement with ex¬ 
periments (^ . Then, using the obtained formulae for 
aij, g, and M, they can be automatically obtained as 
c^xx = ctyy = 14.1, ttxy = Q, g = 58.7, and M = 2.08 
for B II trigonal, which agrees well with the experimen¬ 
tal value of 2.12. The only band that can increase M 
is T^(2) (Fig. 1 in SM). Therefore, and surprisingly, a 
band leV far from the band in which carriers resides can 
enhance the magnitude of M by a factor of two. This 
large interband effect provides a quantitative solution to 
the second puzzle. 

An additional cross-check is provided by the evolution 
of M with alloying or by applying pressure. One 
may naively expect that substituting bismuth with anti¬ 
mony can be assimilated to chemical pressure and there¬ 
fore the results should be identical to applying physical 
pressure. But there are two essential differences. First, 
the substitution changes the strength of the SOI, while 
the pressure does not. Second, the lattice structure be¬ 
comes more rhombohedral by substitutio n 12911 . while it 
tends to approach cubic by the pressure [30|. Figure [3] 
shows the Nernst signal of Bi as a function of EB~^ for 
pure Bi at P = 0, Bio, 96 Sbo.o 4 alloy at P = 0, and pure Bi 
at P = 1.37 GPa. (See SM for details on experiments.) 
We can distinguish the peaks of n = 0"'' and 2~ from 
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FIG. 3. Nernst signal as a function of FB~^ at different 
temperatures for (a) pure Bi at P = 0, (b) Bi with 4% Sb 
substitution at P = 0, and (c) pure Bi at P = 1.37 GPa. 
The magnetic held is along the trigonal axis and P is the 
frequency of the quantum oscillation. Peaks are identihed by 
their Landau sub-level indexes. 


their different temperature dependence. In the case of 
Sb substitution, the 0+ and 2“ peaks shift, but the or¬ 
der is unchanged. Applying pressure, on the other hand, 
inverts the order. This difference can be seen by plot¬ 
ting M as a function of Sb content and pressure shown 
in Figs, m (a) and (b), respectively. M increases with 
substitution but decreases with pressure. Note that in 
both cases as shown in Figs. |4] (c) and (d) the frequency 
of quantum oscillations decreases as it was found before 


31L l32| . This implies that the Fermi surface shrinks be¬ 


cause the overlap between the conduction band at the 
L-point and the valence band at the T-point is reduced. 


These experimental results can be naturally inter¬ 
preted by the present theory in terms of the relevance 
of the interband contributions. Figures |4] (e) and (f) 
show the theoretically obtained energy differences from 
the hole band {T^^) with respect to the Sb content 
and the compression in volume, —Afl/rioj resp ectively. 
(—Afl/flo = 0.041 corresponds to 1.5 GPa [27l|.') In the 


Sb substitution 


Pressure 



Sb content (x) Pressure (GPa) 

o!oO 0.01 0.02 0.03 0.04 0.05 

-AQ / Qo 


FIG. 4. The ratio M, the quantum oscillation frequency, and 
the energy differences from the hole band (T 45 ) as a func¬ 
tion of Sb content (left column) and pressure (right column). 
Experimental and theoretical results are shown by symbols 
and solid lines, respectively. For pressure, the theoretical re¬ 
sults are obtained with respect to the compression in volume, 
—Afl/flo. (—Afl/flo = 0.041 corresponds to 1.5 GPa [ 2 ^). 
Insets of (c) and (d) are the energy shift of the T-point hole 
band, the conduction and valence band at the T-point, taken 
the origin of the energy as the position of the Fermi energy 
of pure Bi at zero pressure. 


case of the substitution, both the higher energy Tg''(2), 
T^(l) and the lower energy Tg''(l) go up, so that the 
contribution from the lower energy band increases re¬ 
sulting in the enhancement of M shown by lines in Fig. 
0 ] (a). In the case of pressure, on the other hand, both 
the higher and lower bands move downward, so that the 
contribution from the higher energy band increases re¬ 
sulting in the decrease of M shown in Fig. 0 ] (b). Even 
though the energy shifts are very small, they make a siz¬ 
able change in M. In both cases, the overlap between the 
T-conduction and T-valence bands decreases as shown in 
the insets of Figs. 0 ] (c) and (d). As seen in the figure, 
there is a satisfactory agreement between theory and ex¬ 
periment. For Sb substitution, we adopt a simple virtual 
crystal approximation [s^. As for pressure, we assume 
that the overlap integrals scale with (i“^, where d is the 
bond length [ 2 ^. (See Sec. V and VI in SM for details.) 
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In summary, we studied both experimentally and theo¬ 
retically the ratio of the Zeeman splitting to the cyclotron 
energy M, which characterizes the effect of crystalline 
SOI. A general theory on M was newly derived based 
on the multiband k ■ p theory. By this multiband ap¬ 
proach, we succeeded in solving the longstanding mystery 
for hole-like carriers at the T-point of bismuth, the large 
and anisotropic M, which have not been explained by the 
previous two-band approach. Our results strongly sug¬ 
gest the surprisingly large impact of interband effect of 
the crystalline SOI. Moreover, we gave new experimen¬ 
tal results on M by substituting antimony and apply¬ 
ing pressure, both were quantitatively explained by the 
present multiband k ■ p theory. The quantitative agree¬ 
ment with experiment is another success for the present 
multiband k ■ p theory. Beyond the specific case of bis¬ 
muth, this general scheme can be applied to all materials 
where SOI plays a relevant role, such as group V semimet¬ 
als (Sb, Ab), IV-VI narrow gap semiconductors (PbTe, 
PbSe, SnTe, etc.), and the topological insulator (Bi 2 Se 3 ) 
in which M ~ 2 has been reported [3, [l3,13 • 
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